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ABSTRACT 

We provide a systematic numerical and analytical study of Klein-Nishina 
(KN) effects in the spectrum produced by a steady state, non-thermal source 
where rapidly accelerated electrons cool by emitting synchrotron radiation and 
Compton upscattering ambient photons produced outside the source. We focus 
on the case where q, the ratio of the ambient radiation field to source magnetic 
field energy densities, significantly exceeds unity. We show that the KN reduction 
in the electron Compton cooling rate causes the steady-state electron spectrum 
to harden at energies 7 > 'Jkn, where 'Jkn = l/4eo and eo = huo/^eC^ is the 
characteristic ambient photon energy. This hardening becomes noticeable in the 
synchrotron radiation from electrons with energies as low as O.I'Jkn and changes 
the synchrotron spectral index relative to its Thomson limit value by as much 
as Aa ~ 0.75 for q ^ 1. The source synchrotron spectrum thus shows a high- 
energy "bump" or excess even though the electron acceleration spectrum has no 
such excess. In contrast, the low-energy Compton gamma-ray spectrum shows 
little distortion because the electron hardening compensates for the KN decline 
in the scattering rate. For sufficiently high electron energies, however, Compton 
cooling becomes so inefficient that synchrotron cooling dominates - an effect 
omitted in most previous studies. The hardening of the electron distribution thus 
stops, leading to a rapid decline in Compton gamma-ray emission, i.e., a strong 
spectral break whose location does not depend on the maximum electron energy. 
This break can limit the importance of Compton gamma-ray pair production on 
ambient photons and implies that a source's synchrotron luminosity may exceed 
its Compton luminosity even though q > 1. We discuss the importance of these 
KN effects in blazars, micro-quasars, and pulsar binaries. 
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1. INTRODUCTION 

Probably because of the initially greater sensitivity of radio telescopes, studies of the 
emission from non-thermal sources, e.g., the powerful radio jets found in AGN, focused on 
synchrotron radiation from relativistic electrons moving in a magnetic field. It was quickly 
realized, however, that the same relativistic electrons would also Compton upscatter any am- 
bient low-energy photons to produce emission at much higher, e.g.. X-ray and gamma-ray, 
energies. A classic discussion of the emission spectrum expected from a gas of relativis- 
tic electrons where synchrotron radiation and Compton scattering are the dominant energy 
loss mechanisms may be found, for example, in Felten & Morrison (1966). Several of the 
standard approximations still used today are presented there, e.g., the delta function ap- 
proximation which provides a one-to-one relation between the energies of the electrons and 
the synchrotron and Compton upscattered photons produced by them. Of particular inter- 
est, those authors note that the processes of Compton upscattering in the Thomson regime 
and synchrotron radiation both lead to quasi-continuous electron energy losses that are pro- 
portional to the square of the electron Lorentz factor. (This is in fact not a surprise given 
that synchrotron radiation may be viewed as the Compton upscattering of virtual magnetic 
field photons [see Blumenthal & Gould 1970].) This implies that the shape of the energy 
distribution of cooled electrons does not depend on which process dominates their cooling 
and that the ratio of the luminosities for the resulting Compton and synchrotron emission 
components simply goes as Lc/Lgyn = Uq/ub. where Uq and Ub are respectively the (co- 
moving) low-energy radiation and magnetic field energy densities inside the source. This is 
very useful to know, for example, since if one measures Lc and Lsyn and knows mq, e.g., 
if the relevant low-energy ambient photons are due to the cosmic microwave background, 
then one can immediately derive the magnetic field strength in the source. Also, the shapes 
of the synchrotron and Compton emission spectra will be similar, with the Compton spec- 
trum simply shifted up in energy relative to the synchrotron one by the factor co/es where 
es = B/Bcr, Bcr = 4:A x 10^^ G, and eo = hvo/irieC^ is the typical ambient photon energy. 
(For mono-energetic ambient photons, the spectral shapes are in fact identical except for the 
effects of self-absorption in the low-energy portion of the synchrotron spectrum.) 

It is important to remember, however, that all of these convenient rules of thumb, used 
in many interpretation papers, break down when the electrons are sufficiently energetic to 
scatter with ambient photons in the Klein-Nishina (KN) limit, i.e., when 7 > 7xAf = l/4eo 
(e.g., Blumenthal & Gould 1970). In this case, the relation between upscattered and incident 
photon energy changes (the upscattered photon cannot have more energy than the incident 
electron), electrons can lose most of their energy in a single scattering rather than cooling 
quasi-continuously as a result of many small energy losses, and to the extent that a cooling 
rate is relevant, the ratio of Compton to synchrotron energy losses is now a decreasing 
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function of energy and is less than uq/ub- A careful re-examination of the implications of 
treating the Compton scattering process correctly in this limit is timely now that satellite 
and ground-based gamma-ray telescopes have shown many extragalactic and some galactic 
compact objects to be strong emitters of GeV and TeV radiation. If produced by Compton 
up-scattering, i.e., the Inverse Compton (IC) process, then these gamma-rays likely result 
from scatterings in the KN regime, and we must be accordingly careful in the interpretation 
of the multi-wavelength observations for these objects. 

Some of the possible KN "corrections" have already been discussed or are well-known. 
An obvious limit is the case when the ambient (target) photon energy density in the source 
is small and the energy density ratio q = Uq/ub is thus ^ 1. In this case, synchrotron losses 
always dominate over Compton losses regardless of electron energy. If we have a source where 
electrons are rapidly accelerated to high energies and subsequently cool via radiative losses, 
then the steady state spectrum of the source (assuming on-going acceleration) is determined 
solely by the magnetic field and the "injection" energy spectrum of electrons produced by the 
rapid acceleration process. We can therefore consider the electron distribution to be fixed, 
e.g., a power law with oc 7"*, and just need to carry out the relatively easy computation 
of the Compton spectrum upscattered by such an electron power law. For simplicity, let 
us assume that the ambient radiation field is approximately mono-energetic, again with 
characteristic energy eo. In this case, the resulting Compton emissivity is approximately a 
broken power law, with the usual Thomson-limit result of je oc e"^^"^"^^"^ for e <^ exN, and 
the extreme KN-limit result of oc e~'^log(e) for e ^ ekn where ekn — Ikn, e.g., see 
Aharonian & Atoyan (1981). In other words, the synchrotron spectrum is unchanged from 
a Thomson-approximation calculation while the Compton spectrum shows a strong spectral 
break at e ~ eKN with a change in spectral index Aa = {s + l)/2. For s > 2, which seems 
typical for the known very high-energy gamma-ray sources, Aa > 1.5, i.e., is large, and 
given the usually poor statistics at the highest energies, this spectral break can easily be 
misinterpreted as a an exponential cutofi^ due to the maximum electron energy in the source. 

A much less obvious but still important limit is the case q ^ 1. For uq sufficiently large, 
one can in principle go the opposite limit and ignore the effects of synchrotron cooling, solv- 
ing directly the kinetic equation for the evolution of electrons due to Compton scattering. 
The steady state electron and upscattered photon spectra obtained for continuous electron 
acceleration and complete electron cooling via Compton scattering are discussed in detail by 
Zdziarski (1989) (see also Zdziarski & Krolik, 1993). The main conclusion is that the equi- 
librium electron distribution hardens for ^ > ^kn compared to the Thomson approximation 
result. This is because electron energy losses are relatively less efficient in the KN limit and 
electrons remain longer at higher energies. This hardening of the electron spectrum compen- 
sates the decreased efficiency of Compton upscattering, and the resulting photon spectrum 
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shows no "Klein- Nishina break" at e ~ eic,KN-- Interestingly, for a rapid acceleration process 
that effectively injects electrons into the source with an energy distribution (5(7) oc 7"^ with 
p = 2, the Compton emissivity is oc e~^, i.e., exactly the answer that would have been 
obtained in the Thomson limit. Note, though, an important difference from the Thomson 
result is that the correctly computed spectrum always cuts off at e ~ 7max, the maximum 
electron Lorentz factor, and that this cutoff is independent of the target photon energy eo- 
Also, for p 7^ 2, the Compton spectrum in the KN regime does not follow the slope of the 
spectrum in the Thomson regime, steepening a bit if p > 2, and hardening if p < 2. 

The analysis just mentioned, however, do not include synchrotron losses and do not 
consider the synchrotron emission produced as a result of these losses. (The condition q ^ 1, 
in fact only guarantees that synchrotron losses are negligible for 7 < Jkn- At higher energies, 
the Compton loss rate decreases due to KN effects and eventually synchrotron coohng always 
dominates.) They thus miss several effects that are potentially important for high energy 
sources. First, the hardening of electron distribution in the KN regime leads to a hardening 
of co-spatially produced synchrotron radiation. Prior work which touches on this aspect 
includes: Dermer & Atoyan (2002) which invokes KN effects to help explain the production 
of X-rays in large scale jets via synchrotron radiation; Ravasio ct al. (2003) which invoke KN 
effects to explain the spectral "glitch;" Kusunose, Takahara & Kato (2004) and Kusunose 
& Takahara (2005) which study KN effects in the context of FSRQ (fiat-spectrum-radio- 
quasars). Secondly, for sufficiently high electron energies, synchrotron losses dominate and 
the electron distribution "saturates" to one with the same slope as in the Thomson regime, 
but with a normalization factor q times larger. This effect was included by Khangulian & 
Aharonian (2005) in their studies of outflows from compact objects in high mass X-ray binary 
systems, but not by Kusunose & Takahara (2005), even though some of their calculations 
involve electron energies high enough for synchrotron losses to be important. 

In this paper, we provide a systematic study of KN effects in steady-state sources with 
g 3> 1, covering all the effects mentioned above. We use accurate numerical techniques 
that solve the exact integro-differential equations for the steady-state photon and electron 
energy distributions. However, we also present approximations that allow one to follow the 
various effects analytically using simple algebraic functions. We focus our studies on the 
case where the ambient radiation field is dominated by external photon sources with mono- 
energetic or power-law spectra. The paper is organized as follows. In §2 we analyze inverse- 
Compton electron energy losses in the KN regime and compare them with the corresponding 
synchrotron energy losses. In §3 we present general and approximate formulas for the inverse- 
Compton emissivities for the cases of isotropic and beamed ambient radiation fields. The 
energy distribution of relativistic electrons in a steady-state source with continuous electron 
injection is discussed in §4. The KN effects in the synchrotron and inverse-Compton spectra 



- 5 - 



produced by such a source are analyzed in §5. Our results are discussed in the context of 
specific astrophysical sources in §6 and summarized in §7. 



2. ELECTRON ENERGY LOSSES 

2.1. Inverse Compton energy losses 

The rate of inverse Compton (IC) energy losses of relativistic, isotropically distributed 
electrons is (see Appendix A) 

n\ic = 7, nUol Fkn , (1) 



3meC^ 



where 

"to 



- j fKN(b)u^od^o : (2) 



uq = p"''"''"^ u^,.deo is the total energy density of the radiation field, u^,. is the energy dis- 
tribution of the ambient photons, and wc rc-cxprcss the electron energy as 6 = A'-fto noting 
6=1 corresponds to the transition between the Thomson and KN scattering regimes. The 
function fxNip) is given by Eqs. (A12) and (AlO) in Appendix A. For 6 <^ 1 (Thomson 
limit), Jkn ~ 1; for 6 > 1 (KN limit), Jkn ^ {^/{2p)){\nh - 11/6). For h < 10^ fKN{h) 
can be approximated by 

fKN ^ • (3) 

(1 + 5)1.5 ^ > 

To estimate the effects of having an extended target photon energy distribution, let us 
assume that is a power law oc eg" and that h^ax — ^^Q,maxlmax < 10^, so we may use 
approximation (3). Then we may write, 

Fkn « / e-"°c^eo + / eo """'-'deo, (4) 

where eq^kn — 1/47. effects become important and Fkn < 1 for 7 > 7xAr = l/(4eo,maa;)- 
For CKo < —0.5, Compton energy losses for electrons with 7 > ^kn are dominated by scat- 
terings on photons with the highest energies, ~ ^o^max- Wc may then treat the photon dis- 
tribution as mono-energetic, with Fkn — fKNip = ^l^.Q^raax)- For —0.5 < ao < 1, Compton 
losses are instead dominated by scattering on photons with energy ~ eq,kn- In this case, we 
may make the so-called "Thomson-edge" or "Klein-Nishina cutoff" approximation that 

Fkn ^ fir: ^-ao^^ ^ (eo,xiv/6o,n^a.)-""+^ = (5) 
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if rewrite the electron energy as 6 = 4eo,maa;7 ^^^d have eo,mm <S ^o,max- A useful approxima- 
tion that interpolates to the Thomson limit for 6 — > and has good accuracy for < ckq < 1 
is then given by, 

If we instead have cto > 1, the Compton losses of the electrons are dominated by scatterings 
on the lowest energy photons available, at ~ eo.mm- We may then use the approximation, 
Fkn — fKN{b — 4:'yeo,min)- If we further have 4:'yrnax^o,min < 1, then in fact all scatterings are 
effectively in the Thomson regime even if higher energy ambient photons are present and 

bmax > 1- To gauge the accuracy of approximations (3) and (6), wc compare in Fig. 1 the 
approximate values of FKN{b) to the exact ones for three different ambient photon spectra: 
mono-energetic, and power laws with «o = —dlnUf^Jdlneo = 0.0 and 0.5. From an analysis 
similar to that for the power law case, one can show that F^n for a Planckian (thermal) 
distribution is well-approximated by treating the Planckian as a mono-energetic photon 
distribution with energy 2.8kT. 



2.2. Synchrotron energy losses vs. inverse-Compton energy losses 

Noting that the rate of synchrotron losses in a tangled magnetic field of strength B is 

n\syn = 2^Bl 

where ub — B^/Stt is the magnetic energy density, we have (see Eq.l) 

7/c 



Isyn 



qFkn (8) 



where q = uq/ub- Since Fkn < 1 for any value of b, for g < 1 the energy losses for all 
electrons are dominated by synchrotron radiation. For q > 1, the energy losses for electrons 
with 7 < 75 are dominated by inverse Compton scattering while losses for electrons with 
7 > 7s are instead dominated by synchrotron radiation, where 7^ = bs/4:eo,max, and bg is the 
solution of the equation qF^N — 1, plotted in Fig. 2. For a mono-energetic ambient photon 
spectrum (see Eq.3) 

bs ~ g'/' - 1 , (9) 
while for a power-law spectrum with < ccq < 1 (see Eq.6), 

bs ~ - 1 . (10) 
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The relative role of the IC and synchrotron energy losses is illustrated in Fig. 3. We 
show there 

iic qFKN{h) 



jtot 1 + q^KNib) ' 

and 



jtot 1 + qFKN{b) ' 

for three different spectra of the ambient radiation field, assuming that jtot — iic + isyn- 



(11) 

(12) 



3. INVERSE-COMPTON EMISSIVITY 

The inverse Compton emissivity for isotropically distributed electrons is 

fdP^IC)h) , /ION 

Je{ic) = / n^dj, (13) 



where 

&P,,a,(l) _ SNUe,-,)^^^^, ^ (14) 



dn dedQ 

is the IC-power per unit photon energy per solid angle per electron. For two particular cases, 
mono-directional beams of target photons and an isotropic radiation field, the respective 
scattering rates are given by Eqs. (Al) and (A4). Inserting them into Eq. (13) gives, for the 
beamed mono-energetic external radiation field, 

eie(7C)(^) = Y^ccttMo (^^^ J ^f{j,e,eo,e)dj , (15) 
and for the isotropic ambient radiation field, 

eje(7C) = Y^carMo (^^^ j^fiso{l,e,eo)d'y, (16) 

where 9 is the scattering angle (the angle between the photon beam and the direction to the 
observer), the functions /(7, e, eo, 9) and fiso{j, e, eo) are given by Eqs. (A3) and (A7). These 
formulas can easily be generalized for any spectrum of external radiation field by replacing 
Uq by Meofieo and integrating them over cq. 

A useful first approximation for the IC spectrum (e.g., see Coppi & Blandford 1990) 
may be obtained by making the delta function approximation /(e,7, eo) oc S[e — e(/c)(7)]) 
where 

_ Jef{-f,e,eo,9)de 
J fu, e, eo, 9)de 
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and 

J fisoKl, e, eo)de 

are the average energies of photons produced by the scattering of photons of energy Cq by 
electrons with energy 7. A useful quantity to know when dealing with the KN regime is the 
electron's scattering inelasticity, i.e., the average fraction of the incident electron's energy 
that is transferred to the scattered photon: A — e/c'(7)/7- This is shown in Fig. 4 for the 
case of a mono-energetic ambient photon field. 

Inserting the delta function approximation into Eqs. (15), and (16), we can rewrite the 
IC emissivities as 

«,,c,W.!!22|,|,,,„.e^x(7.^)5^^. (19) 

and 

, . , 9 din 7 _ 

Qe{ic) - -j^ \l\icmec' -p^ , 20 

VK alne(/c) 

respectively, where 

j /iso(e,eo,7)eae e/c(7) 

Note that in the Thomson limit, lie ^ 7^, so that dln^/dlnejc — 1/2, while in the KN 
limit, ejc ~ 7 so that dln^/dlne(^ic) — 1. 

We can make an analogous approximation for the synchrotron emissivity, provided that 
the magnetic fields are isotropic and 7^00: < Bcr/{^B). We have then, 

n^7 , . I 2 c?ln7 1 n^7 2 (r,r,\ 

^Hsyn) - \l\synmeC ^ = 2 17 \l\synmeC (22) 

where Isyn = (4/3)7^(5/5^^) and = 27rm^c^/(/ie) ~ 4.4 x lO^^Gauss. 



4. STEADY-STATE ELECTRON ENERGY DISTRIBUTION 

Under the assumptions that the time scales for the acceleration of individual electrons 
are much shorter than time scales of their energy losses, that electrons do not escape from the 
coohng region, and that they also do not suffer adiabatic losses, the evolution of the electron 
energy distribution, n^, can be described by the integro-differential equation (Blumenthal & 
Gould 1970) 

^ = -^(n,|7l)-n,y^ C{^,i)di + j n,C{i,^)di + Q, (23) 
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where I7I is the energy loss rate due to loss processes that may be approximated as being 
continuous, e.g., synchrotron radiation, C(7, 7') is the probability per unit time for Compton 
scattering of an electron with energy 7 to energy 7', and Q is the electron injection function, 
i.e. the electron production rate per unit time, per energy and per volume. The transition 
rates (7(7,7') have been derived by Jones (1968) for a mono-energetic ambient radiation 
field, and by Zdziarski (1988) for power-law ambient radiation fields. 

To illustrate in a simple way the effects of KN corrections, we will use a continuity 
version of the kinetic equation 

^--^Kl7U)+g, (24) 

despite the fact that in the KN regime the fractional electron energy losses per scattering 
are not negligible (see Fig. 4). Such a simplification can be justified by noting that the 
results obtained by Zdziarski (1989) show that the electron distributions obtained by using 
the exact integro-differential equation versus the continuity equation are qualitatively very 
similar and that significant differences (up to a factor few) occur only if both the electron 
injection function and the ambient radiation field are mono-energetic. In the cases studied 
here, the differences are further reduced by the increasingly important role of the (continuous) 
synchrotron energy losses at the highest energies, particularly if hmax > bs- We demonstrate 
this explicitly in Fig. 5-7 where we compare the results of exact calculations made using 
the code of Coppi (1992) with the approximate results obtained by treating Compton losses 
continuously in that code. For many applications that do not involve fitting data and thus 
do not require very high accuracy, the continuous energy loss approximation should suffice. 

A steady-state solution of the continuity equation, dn-y/dt — 0, is 

where we assume that electron energy losses are dominated by inverse Compton process and 
by synchrotron radiation, i.e. that 

4:CaTUB 2/ 7-> \ 

|7tot| = Hsynl + nic\ = ^2 ^ (1 + qFkn) ■ (26) 

For bmax ^ bs and q ^ 1, the steady-state then distribution has two asymptotes: one for 
6 <IC 1, where electron energy losses are dominated by Thomson scatterings, i.e. F^n = 1; 
and one for bs <^ b <^ bmax, where electron energy losses are dominated by synchrotron 
radiation, i.e. qFkn <^ 1- The respective equilibrium electron distributions are: 7^n-y oc 
J Qd^/{1 + q) and 7^% oc J Qd^. For 1 < h < bg, despite the decrease of the KN 
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cross section, electron energy losses are still dominated by the IC process and the electron 
distribution is 7^n-y oc j^Qd'y / {qFKN}- 

For a rapid acceleration mechanism that produces a power law electron injection spec- 
trum, Q oc 7~^, with p > 1, the two asymptotes of the steady state distribution (for 6 <^ 1, 
and h ^ hs) are power laws with the same spectral index, s — p+l, where s = —dlnn^/dln.^. 
For 1 <^ 6 <^ 6s, i.e., for ^kn < 1 <1s-i the electron distribution is /iarderthan in the asymp- 
totic regions, with a spectral index reaching s = p+ 1 + As, where As ~ dXnFKN / d^^^l < 0. 
Using the analytical approximations for Fkn given in §2 and going to the limit q ^ 1, one 
finds that for ambient radiation fields with mono-cncrgctic or very sharply peaked photon 
energy distributions. As ~ —1.5, while for softer external fields, with ao > 0.0, As ~ ao ~ 1- 
When q < 10^, the maximum hardening. As, of the electron distribution is a function of 
q and decreases with decreasing q. (See Fig. 8 for an example of how the corresponding 
hardening of the electron synchrotron depends on q.) Note that to the extent the continuous 
Compton cooling approximation is a good one, this result is independent of the injection 
index p. Of course, for bmax < bs the electron energy distribution will not have the maximum 
deviation possible from the Thomson-limit (6 <^ 1) asymptote and As may not reach its 
maximum value. In this case, the KN induced bump or "excess" at the high energy end of 
the electron distribution is correspondingly less prominent. This is demonstrated in Fig. 5, 
where the electron distribution is computed for four different values of bmax- 



5. SPECTRA 

5.1. Klein-Nishina effects 

Examples of the steady-state electromagnetic spectra produced for the case of power- 
law electron injection function are shown in Figs. 6 and 7. The IC spectra are computed 
assuming an isotropic ambient radiation field. The characteristic photon energies marked on 
the figures have the following definitions (see Eq. 18): eic,KN = ^ic{ikn)] ^ic,s = ^icils)', 
^syn,KN = ^syn{lKN)\ (7s); and e+ = 2/eo,maa;, which is the approximate threshold 

energy for photon-photon pair production on the ambient photon distribution. 

As one can see, the high energy portions of the inverse-Compton (IC) and synchrotron 
spectra behave very differently. The IC spectra do not change very much after crossing 
^ic,KN- This is because for g ^ 1, Compton scattering is still dominating the cooling of 
the electrons responsible for IC photons at e > ejcKN- Even though those electrons scatter 
in the KN regime, the decreased efficiency of Compton scattering is compensated by the 
corresponding increase in their equilibrium density, as discussed in Zdziarski & Krolik (1993). 
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As one moves to higher photon energies, however, synchrotron coohng becomes relatively 
more important for the electrons responsible these IC photons. Eventually, when the energy 
of the relevant electrons exceeds 7^, synchrotron cooling rapidly dominates and the fraction 
of the electron's energy going into the Compton component plummets. The result is a sharp 
steepening or break of the IC spectrum at energy eic,s, which is independent of the maximum 
electron energy '^max > Is- 

When one instead considers synchrotron emission, the hardening of the electron dis- 
tribution at 7 > jkn is directly reflected in the synchrotron spectrum. The result is a 
synchrotron spectrum that can harden dramatically at e > esyn.KN, forming a synchrotron 
"bump" if bmax < bs- (Note that if g ^ 1, the hardening of the synchrotron component 
can already be noticeable at even lower energies, ~ ^syni^-^lKN)-) For b^ax ^ bg, the syn- 
chrotron spectrum at e > esyn,s asymptotes to a spectrum with the same slope as the low 
energy (Thomson regime) asymptote but with a normalization that is a factor q higher. 

An interesting consequence of the very different behaviors of the high energy portions of 
IC and synchrotron spectra is that for the case of a flat electron injection spectrum {p < 2) 
with b^ax > bs, the luminosity of the synchrotron peak - located around esyn,max = ^synilmax) 
- will exceed the luminosity of the IC peak - located around eic,s, no matter how large we 
make q. This dominance of the synchrotron component, even though q — 30 ^ 1, is 
demonstrated in Fig. 7. 

All the spectral features just described can be reproduced analytically, by using the 
equations 



which are obtained after insertion of Eq. (25) into Eqs. (20) and (22), respectively. For 
Imax ^ 7s and q ^ 1, the latter implying bs ^ 1 i.e. 7^ ^ "yxN, the IC and synchrotron 
spectra can be characterized as the superposition of three components produced by electrons 
with: 7 < ^kn; Ikn < 7 < 7s; and 7^ < 7 < ^max- These components are: 




(27) 



and 




(28) 




(29) 
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where 7 ~ •>/(3/4)(e/eo,r„aa;) if e < e/cftriv, and e ~ 7 if e > e/ci^rjv; and 

4^ . , , /■ f V(2g) if e < e.j/n,/fiv 

j.(.yn) J Qd^^ ) i/{2qFKN) if e,j,„,,^iv < e < e,j,„,, , (30) 



I 1/2 if esj/n,s < e < e 



syn,max 



where 7 = v/3eficr/(4fi). 

For electron injection Q oc 7"^ and p > 1, the various components are power laws, i.e., 
we have: 

4^ . r if e < e/cifiv 

^ oc \ e-(f-i) if ejc/fiv « e « e^c,. , (31) 

U-(^-^+^«) if e,c,. « e « e,c,n.a. 
where 5q; ~ 1.5 for ckq < —1, and 5q; ~ 1 — ckq for < ckq < 1, while 



4vreje(sj/n) 

oc 



e if e < e^j,^,/^ jv 

,-W2+Aa) if e,j,„,;f jv < e < e,j,„,, , (32) 



e-P/^ if e,„„,5 < e < e. 



~syn,s >^ ^ '^syn,max 

where Aa = — (dln(l/Fft:iv)/c^lne)|5(xei/2- For cto < —1, Act ~ —0.75, while for < ao < 
1, Aa ~ — (1 — q;o)/2. We should emphasize here, that estimation of the hardening of 
synchrotron spectrum slope in the band [esyn,KN') ^syn,s\ is very crude and corresponds to its 
maximum value, which as noted previously, can be reached only for q > 10'^. As shown in 
Fig. 8, the synchrotron spectral hardening will be weaker for lower values of q. 

Although we have only presented and discussed emission spectra computed for an 
isotropic ambient radiation field, we would like to emphasize, that with our generic as- 
sumption of isotropic distribution of electrons, the electron energy losses, as well, as the 
distribution of electrons (Eq. 25) and their synchrotron emissivity (Eq. 22), do not depend 
on angular distribution of ambient radiation field. Hence, in the case of a beamed ambient 
radiation field the IC spectrum can be simply computed from Eq. (15) or (19) by using 
there the electron distribution given by Eq. (25). In Fig. 10 we show such spectra for three 
different scattering angles and compare them with the one computed for the isotropic am- 
bient radiation field. One can see that the deeper one is in the KN regime, the weaker is 
anisotropy of the scattered radiation. Suppression of the anisotropy is caused by the recoil 
effect. 



5.2. Further effects that may modify the observed spectrum 



To highlight the spectral effects caused by the modification of the electron energy dis- 
tribution due to Compton scatterings in the KN regime, we have considered only spectra 
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produced in the fast electron cooling regime, and we have ignored possibly important pro- 
cesses such as photon-photon pair production and synchrotron self-Compton radiation. We 
also did not consider effects due to the relativistic propagation of the source, which is im- 
portant in objects like blazars. We discuss below how our results can be affected by the 
inclusion of some of these complications. 

The fast vs. slow cooling regime 

Assuming that electron energy losses are dominated by radiative processes, the cooling 
time scale for an electron of energy 7 is (see Eq. 26) 

I7I AarUB'jil + qFkn) ' 

For sources with a finite (comoving frame) lifetime tq, the only electrons that have time to 
cool significantly are those with < tq. Looking first at the Thomson limit (7 < ^kn) of 
this expression, one finds the usual result that only electrons with 7 > 7c have time to cool, 
where 

AuTUotq 

for q ^ 1. The number density of electrons at a given energy that can accumulate in the 
source is roughly (^(7) x niiii[tc(7), ^q]- For power law electron injection, the electron energy 
distribution of electrons therefore hardens by As = 1 for 7 < 7c, leading to a hardening of 
the synchrotron or Compton spectrum by Aa = 0.5. 

For 7 > 7KAr, ? 1, and a mono-energetic ambient photon distribution, the the cooling 
time increases for 7 > '"yxN-, reaching a local maximum at 7 f« 7^, when synchrotron cooling 
begins to dominate. Efficient KN cooling therefore requires tc(7s) < iq- From Eq. (33), we 
have 

3meC 3meCeo q 

Hence, fast cooling in the KN regime requires 

f < (36) 
or, equivalently using bs ~ q'^^^ for 1 -C g < 10^, 

,<(?^)\ (37) 
Of course, if we have b^ax < bg, the upper limit on q is correspondingly weaker. 
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For a power law photon distribution with < cto < 1, we may use Eq. (33) and 
approximation (6) to show that tc ~ 7~"°. For ao > 0, this is a monotonically decreasing 
function of 7, so fast cooling throughout the KN regime is simply guaranteed by the condition 
tcilKN) < tq, or equivalcntly jkn > 7c- This can be translated into the following requirement 
on the energy density of the external radiation field, 



For a photon distribution that is harder than = oi' softer than = 1, the condition for 
a mono-energetic radiation field, i.e., we require tc(7s) < tg- 

Photon-photon pair production 

Because (a) the cross-section for photon-photon pair production is similar in magnitude 
to that for Compton scattering, and (b) the photon threshold energy for pair production, 
e+ ~ i/^o,max, is almost the same as the electron energy ■Jkn ^ ^/^o,max, it is often stated 
that strong pair production is unavoidable in sources where KN effects are important. This 
is not always true, however. 

First, in many applications a better estimate for the threshold energy is in fact e+ = 
2/eo,maa;- Morcovcr, KN effects arc actually important at energies well below l/eo,max, i-e., at 
7 < iKN = l/4eo,max- Furthermore, from Fig. 4., A{jkn) ~ 0.1 — 0.3, so that an electron of 
energy jkn Compton actually upscatters photons to typical energies eic,KN = -^knIkn 
e+. Taking A at larger energies to be ~ 0.5, we see that we in fact need a source with 
bmax > = l+llKN ~— 16. Siuce KN distortions of the electron spectrum already produce 
visible distortions in the synchrotron spectrum for electron energies h ~ 0.1, this means 
there is a factor ~ 100 in electron energy for which KN corrections are important but pair 
production is not possible. Since Csyn oc b^, this corresponds to a factor 10^(!) in synchrotron 
frequency. 

Second, even if we have bmax ^ b+, pair production may still not be important. When 
one includes the effects of synchrotron cooling, left out of KN pair production studies such 
as Zdziarski (1988), we have seen that one obtains a strong break at eic,s corresponding to 
the electron energy bg where synchrotron cooling starts to dominate. Thus, if 6+ < bmax but 
b^ > bg, pair production occurs but the luminosity of the pairs that are produced (and the 
spectral distortions they induce) will not be bolomctrically important. For a mono-energetic 
photon distribution, bg — q^^^, and one thus has bg < b+ for any q < 60, independent of 
the actual value of the maximum electron energy, bmax, e.g., see Fig. 6a and Fig. 7. Note, 
though, that for an ambient photon distribution that is not mono-energetic, e.g., a power 
law with ao > —0.5, KN effects are not as strong as for the mono-energetic case because 
lower energy photons that scatter in the Thomson regime are available. The spectral break 



SrUeC 



3meCeo 



(38) 



Uq > 



^iKNCrrtQ 
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due to e/cs therefore occurs at higher energy and pair production can be important for much 
lower g, e.g., see Fig. 6b. 

Even if pair production is energetically possible, the preceding discussion says nothing 
about whether the optical depth to pair production r^^ actually exceeds unity, and it ignores 
the effects of photon anisotropy (which raise the pair production threshold). In a realistic 
source, the extent and geometry of the external radiation field must be taken into as well 
as its absolute intensity. In galactic pulsar wind applications, for example, the ambient 
radiation field due to the companion star is often highly anisotropic in the source region of 
interest. In the most general must be treated as a free parameter. In particular, if 

the source region has an effectively infinite lifetime and is not expanding, which would induce 
adiabatic losses, we are essentially free to choose as low an external field energy density as 
we like without violating the fast KN coohng constraints discussed above. (If we choose 
too low an external energy density, of course, higher order like bremsstrahlung start to be 
important.) If the source is not static, though, there are interesting limits we can place 
on r-y-, if we demand efficient KN cooling. For example, consider a source of size R with a 
characteristic lifetime or expansion timescale ~ R/c. For simplicity, assume also that the 
radiation field is mono-energetic. Then, taking r^^ '2:i Q^UqUtR where no = Uq/cq and using 
the fast-cooling condition of eqn. (35), we have 

r,, > (39) 

This implies, for example, that we can be in the fast KN regime and still have r^^ < 1 for 
g<37. 

Synchrotron- Self Compton (SSC) Effects 

In addition to the external ambient photons, synchrotron photons from the cooling elec- 
trons will always be present in the source. The synchrotron photons spectrum is typically 
broad and extends to low energies, i.e., it is not well-approximated by a mono-energetic 
photon distribution. The SSC spectrum produced by synchrotron photon upscattering can 
thus be quite different in shape from the spectrum produced by the upscattering of external 
photons, in particular it may extend to higher energies. This may be important in certain 
applications. For a very weak external radiation field or a very compact source, the Compton 
losses due the synchrotron photons may in fact dominate over the losses due to the external 
field, leading to significant changes in the equilibrium electron distribution and thus the 
emergent spectrum. These will be explored further in a subsequent paper. The analysis 
is more complicated than in the present case because the target radiation field and the 
electron distribution must be determined simultaneously and self-consistently. (As we have 
seen, many quantities depend sensitively on the details of the low-energy photon spectrum.) 
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Nonetheless, the main effects described fiere still occur, e.g., the synchrotron spectrum is 
harder than expected from a Thomson- limit analysis and for sufficiently high electron injec- 
tion energies, there will be an energy 7^ above which synchrotron losses dominate, leading 
to a corresponding break in the Compton spectrum. 

Relativistic source propagation effects 

If a source is moving with relativistic speed (bulk Lorentz factor F ^ 1), and the 
external radiation field in such a frame is isotropic, and with an energy density Uq peaked 
at photon energies ~ eo, then in the source rest frame, the energy density of the external 
radiation is times larger and the energies of seed photons are F times higher. Then the rest 
frame quantities relevant to our discussion of KN effects are b = 4F7eo, and, in particular, 

— l/(4r7eo), and q ~ T^uq/ub- Since for F ^ 1 the head-on approximation applies, 
the IC Compton emissivity can be computed using Eq.(15) with 

cos^^ = -cos^„,, = -^— (40) 

where ip'^^^^ is the angle between the jet axis and the direction to the observer in the co-moving 
frame while ipobs is the same angle but as measured in the lab frame. The observed source 
flux is then 

e..^., ^ ^^^if^ (41) 

where eobs = e^^/(l + z), V = 1/[F(1 — Pcosipobs)] is the Doppler factor, (11 is the distance 
(luminosity distance for cosmological objects), z is the redshift, and V is the volume of the 
source as measured in its rest frame. The hfetime of the source measured in the lab frame 

is rtQ. 



6. APPLICATIONS 
6.1. Blazars 

The clearest and probably most numerous examples of high-gr (radiation dominated) 
non-thermal sources are the powerful blazars where we think we are seeing emission from 
a relativistic jet oriented towards us. In many of these objects, the bolometric luminosity 
is strongly dominated by 7-rays (von Montigny et al. 1995; Mukherjee et al. 1997). A 
day-week variability time scales suggest that these sources are located at (sub-)parsec dis- 
tances from the center. There, the external radiation field, as viewed in the jet co-moving 
frame, is dominated by the powerful Broad Emission Line (BEL) region radiation. For 



-17- 



uo ~ LBEh/i^T^rBELc) and ub ^ LB/(7r(r6lj)^cP), 

q{r ~ rBEi) ~ ~ 25 '—-y — — — — , (42) 

Ub 4ivs^45 

where is the magnetic energy flux carried by the jet and 9j = R/r is the half-opening 
angle of a jet. 

Let us determine now the values of ^kn and 7s and of the corresponding IC and 
synchrotron photon energies for such an external radiation field, assuming that q ^ 1 
and 7„iax- > 7s- Noting that the energies of broad emission lines peak around ~ lOeV 
(eo ~ 2 X 10~^) and that they are seen in the jet co-moving frame as boosted by a factor 
~ r, i.e. b = 4eo7r, we have 

" 4^ ^ io'(r/io)-^ (43) 

and 

7. = bs^KN ~ 10^(g/30)2/3(r/10)-^ (44) 
These electrons Comptonize external photons up to average energies 

efcViv - Akn^knV ~ 2 x 10%AKN/0.U){V/r) (~ IGeV ...), (45) 

and 

efS,s - AslsT) ~ 6 X 10^(^/0.5) (g/30)2/3(P/r) (~ 30GeV...), (46) 
and produce synchrotron photons with average energies 



^fyn,KN = \^KN^V ^ 2 X 10"^ ^^'""^^^L 3 X lO^^Hz...) , (47) 

= ^'C^iv ^ <l'^'eil,KN - 2 X 10-^(g/30)^/^ /^'^^^^/7m (^ ^ x lO^^Hz...) . (48) 

Hence, sources with g ^ 1, 1 < hmax < and power- law electron injection should produce 
synchrotron "bumps" peaking in the IR-UV spectral band, with the closer a given hmax is 
to 6s, the more prominent the bump. As discussed, the presence of "excess' synchrotron 
emission (the bump) in this case is simply due to KN effects and should not be interpreted 
as indicating the presence of a new electron acceleration component or a hardening of the 
low energy injection spectrum. 



We have just estimated the electron injection parameters that would put a luminous 
blazar into the KN regime studied here. Are there any blazars that actually populate this 
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region of parameter space? Since for even the brightest blazars, EGRET could not detect 
gamma-rays much beyond 1 GeV, the extension of the electron energy distribution into KN 
regime cannot be established using EGRET data. However, quite significant constraints on 
the high energy tails of the electron energy distribution are provided by observations of their 
synchrotron spectra. In the study of Padovani et al. 2003, about half of the objects in their 
sample of powerful blazar objects have synchrotron spectra peaking at z/ > 3 x 10^'^Hz, which 
implies 6 > 1. In most cases the spectra steepen in the UV band, indicating an injection 
function with a cutoff or break at 6 ~ a few, i.e., KN effects may be moderately important. 
There are few of these blazars, though, that have a synchrotron peak clearly located in the 
UV to soft X-ray band and may have b^ax ^ 1- For these objects, assuming the external BEL 
radiation field is dominant with the characteristics described above, we then predict a high 
energy spectral break of the IC component at ejc^s ~ 30 GeV, independent of the observed 
synchrotron cutoff energy (i.e., bmax) provided that it is sufficiently large. If the electron 
injection spectrum for these objects is not unusually soft, i.e., we have p < 3, then as shown 
in Fig. 6-8, for example, a significant of their bolometric luminosity actually emerges in the 
synchrotron component. In particular, their IC luminosity should be comparable to their 
synchrotron luminosity, i.e., they would no^be gamma-ray loud objects (with Ljc/ Lgyn ^ 1). 

It would be very convenient if powerful blazars had gamma-ray spectra extending to 
TeV energies, e.g., they would provide very bright sources that could used to constrain 
the intensity of the extragalactic background light via the absorption of their gamma-rays 
(Coppi & Aharonian 1999). Unfortunately, for the BEL parameters we have used, this is 
impossible unless q is extremely (and implausibly) large. One caveat to this conclusion is 
that we have not included SSC effects in our estimates, and because the synchroton emission 
is much broader in energy than the BEL one, 7^ and ejc^s could move to higher values. Note, 
though, that £+, the energy above which pair production on the BEL becomes possible, is 
only ~ 50 GcV. Strong Compton TcV emission therefore seems unlikely unless it occurs far 
from the BEL and the typical ambient photon energy is in the near-infrared range. 



6.2. Micro-blazars? 

Only two or three EGRET sources have been identified with micro-quasars (Parcdes 
ct al. 2000; Massi et al. 2004; Combi et al. 2004). The fact that these objects arc High 
Mass X-Ray Binary (HMXB) systems containing massive and very luminous companion 
stars strongly suggests a Compton origin for their 7-rays. As in blazars, the EGRET ob- 
servations unfortunately do not provide constraints on ^max for these objects. Furthermore, 
because these sources are relatively weak and completely dominated in the optical/UV band 
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by radiation from the companion stars, their synchrotron component cannot be identified. 
However, GRS 1915+105 proves that XRB have the abihty to produce much more powerful 
and relativistic jets (Fender & Belloni 2004) than in the sources just mentioned. If such a 
powerful jet were to occur in an HMXB and it pointed toward us, we would see a micro- 
blazar, with relativistically boosted non-thermal radiation dominating over all spectral bands 
(Georganopoulos, Aharonian, & Kirk 2002). Since the spectra of HMXB companion stars 
peak at ~ 10 eV, the same value for the BEL in quasars, the values of ^kn and 7^^, and of 
^ic,KN and eic,s are of the same order as for blazars and, therefore, like blazars, they are 
expected to be GeV emitters and not strong TeV emitters. Due to the much stronger mag- 
netic fields in micro-quasar jets, the synchrotron spectral bumps, produced if bmax ^ 1; are 
predicted to peak in the UV/soft X-ray band. Hence, if some of the ULX (Ultra-Luminous 
X-ray) sources are in fact micro-blazars pointed at us, they should be strong 7-ray emitters 
(Georganopoulos et al. 2002). 

However, we must remember, that very large q is available only on size scales comparable 
to those of the binary system. For electron acceleration occurring well down the jet, outside 
the binary system, the companion star radiation chases the relativistically moving emission 
region from behind and its energy density is thus Dopplcr deboosted when viewed in the 
jet frame. (The energy density of the companion star radiation field will be further reduced 
by the usual factor, oc 1/r^ where r is the distance from the binary, but this effect can be 
canceled out by the fact that the magnetic field energy in a conically expanding jet also 
drops as 1/r^.) Closer in to the central object, jets typically have much stronger magnetic 
fields and, therefore, synchrotron radiation will dominate electron energy losses, even in the 
Thomson regime. In this case the IC spectrum is expected to break at ejc\KN, i-e., ~ 1 GeV. 
This would be the case if, as in blazars, the jet energy is dissipated at 10'^^'^(r/rg, ) where Vg 
is the gravitational radius of the compact object, i.e., on scales 10^ times smaller than the 
typical size of an XRB system. 

It should be noted that because the radiation field of the companion star is not symmet- 
ric about the jet axis, detailed computations of the non-thermal spectra from XRB require 
the integration of emissivities given by Eq. (15) over the energy distribution of the exter- 
nal radiation field, taking into account that in the source comoving frame, the radiation is 
boosted by a Doppler factor that depends on the direction of the incoming photons (Khangu- 
han & Aharonian 2005). In other words, the exact geometry of the system is important. 
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6.3. Pulsars in HMXB 

PSR B1259-63 and PSR J0045-73 provide examples of non-accreting binary pulsars 
with massive companions (Johnston et al. 1992; Kaspi et al. 1994); and PSR B1259-63 was 
recently identified as a TeV source (Aharonian et al. 2005). Ball and Kirk (2000) envisaged 
two scenarios for production of high energy radiation in such systems, Comptonization of the 
radiation field of the companion luminous star directly by the pulsar ultra-relativistic wind 
(i.e., a bulk-Compton scenario) or by particles accelerated in the terminal wind shock formed 
due to the confinement of the pulsar wind by the wind from the companion star. Since for 
finite values of g ^ 1, the IC spectrum has a high energy break at eic,br = ^^t^{^ic,s, ^ic,max), 
the condition for efficient TeV production is eic,s > 10® and eic,max > 10®. Since eic,s ~ 
IsbslKN) where ^kn = l/(4eo) ~ 10'^, this translates into the condition bs > 100, i.e., for a 
mono-energetic field, we again need a very large q > 10'^. In the shocked wind scenario, this 
condition provides an upper limit on the strength of the magnetic field in the shocked plasma. 
In the bulk-Compton scenario, the condition is satisfied even for strongly magnetized winds. 
This is because electrons are cold and frozen to the magnetic field fines, hence their energy 
losses via the synchrotron mechanism are negligible and, effectively, the value of q is infinite. 
Of course, to reach TeV energies in this scenario, the bulk Lorentz factor of the wind is 
required to be of the order 10®, which is consistent with estimations of the wind speed in the 
Crab Nebula (Rees & Gunn 1974). 



6.4. Kiloparsec scale jets 

Jets in quasars encounter a variety of radiation fields as they make their way out from 

the central black hole, starting from the radiation field of the black hole accretion disk 
and ending with the cosmic microwave background (CMB). For jets that are relativistic 
from the start, the largest q would be reached near the base of the jet. The lack of bulk 
Compton features in the soft X-ray band suggests this is not the case (Moderski et al. 
2004), and thus that the energetically dominant Compton interactions of a jet with external 
radiation field seem to take place not earlier than in the BEL region, i.e. around 0.1 — 1 
parsec from the center. Large values of q are also expected in sources triggered at 1-10 
parsec distances where the diffuse component of the external field is likely dominated by the 
thermal emission from hot dust (Sikora et al. 2002). At progressively larger distances, the 
jet undergoes Compton interactions with narrow emission lines, with stellar radiation, and 
finally the CMB. As of now there is no direct evidence that at such distances electrons are 
injected with bmax ^ 1, but this may simply be due to the sensitivity and angular resolution 
of present 7-ray detectors. Indirect evidence in favor of bmax ^ 1 and g 3> 1 for kiloparsec 
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scale jets may come from the work of Dermer and Atoyan (2002). They argued that a purely 
synchrotron emission model, with a jet magnetic energy density less than the CMB energy 
density, can successfully explain not just the observed radio-to-optical jet radiation, but 
also the X-ray flux detected by Chandra, despite the fact that the Chandra flux often lies 
above an extrapolation of the radio-optical spectrum. This is achieved by the formation of 
a synchrotron bump above esyn,KN due to the KN effects we have described here. We note, 
however, this particular model requires extremely large jmax, at least one order of magnitude 
larger than -fKN = l/(4recivfB) ^ 2 x 10V(r/10)(l + z), where ecMsiz = 0) ~ lO'^. 

7. CONCLUSIONS 

We have studied, both analytically and using accurate numerical codes, the electromag- 
netic spectra produced by relativistic electrons in a magnetized non-thermal source that is 
immersed in a dense radiation field originating outside the source. We consider the case when 
the steady-state electron energy distribution is determined by the injection energy spectrum 
of the accelerated electrons and their energy losses, dominated by synchrotron radiation 
and Compton scattering which may extend deeply into the KN regime. We concentrate 
on the poorly studied region of parameter space in which the energy density of the radia- 
tion field inside the source exceeds that of the magnetic field, i.e., values of the parameter 
q = uq/ub > 1- Fig. 10 summarizes the three main regions we find for the overall parameter 
space for our problem: 

- In zone I (g < 1), the electron distribution is determined by synchrotron cooling and KN 
effects do not appear in the synchrotron spectrum. The Compton spectrum, however, shows 
a strong break at ejc,KN due to the strong KN reduction in the Compton scattering rate that 
starts for electron energies 7 > '-/kn- Assuming the maximum electron acceleration energy, 
Imax: is Sufficiently large that 6jc,max > eic,KN, the position of this break is independent of 
^ymax- In zone lb (1 < g' < 3), we start to see a hardening of the electron distribution that 
is reflected in the synchrotron emission spectrum. The hardening occurs because Compton 
cooling is now an important contribution to the total electron cooling rate, and Compton 
losses decrease at high energies due to the KN effect. The effect is not large, however. The 
Compton spectrum docs not show significant differences because synchrotron losses start to 
dominate again at 7^ < ■Jkn,- 

- In zones II and III, q ^ 1, and the distortion in the electron distribution due to the 
reduction in Compton cooling is very large. The synchrotron spectrum correspondingly 
shows a strong, hard excess over the Thomson limit asymptote. The Compton spectrum 
below eic,s, however, does not show a strong deviation from the low energy Thomson limit 
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because the KN decline in the Compton rate is compensated by the corresponding increase 
in the electron density. The distinction between zone II and III is that in zone II, the 
maximum electron energy is such that 7max > 7s- Above 7^, synchrotron cooling dominates. 
The energy distributions of the electrons and the synchrotron radiation reach asymptotes 
with the same slopes as in the Thomson regime but with amplitudes enhanced by a factor 
q. The Compton spectrum, on the other hand, shows a strong break at eic,s because the KN 
reduction in the scattering is no longer compensated by a hardening of the electron energy 
distribution. In zone II, then, the Compton spectrum breaks at an energy independent 
of Imax- In zone III, instead, the synchrotron losses never dominate, and the location of 
Compton high energy break is determined by the maximum electron energy just as it is for 
the synchrotron component. The combination of the high energy break and the hardening 
of the synchrotron spectrum at lower energies leads to the formation of a strong synchrotron 
bump. 

The specific conclusions of our study are: 

• The IC spectra have high energy breaks at eic,max, if b^ax < bg, or at e/c.s, if b^ax > bg. 
The former is related to the high energy cut-off of electron injection function, the latter to 
the strong steepening of the IC spectrum, caused by domination of energy losses of electrons 
with 7 > 7s by synchrotron mechanism; 

• Synchrotron spectra undergo strong hardening at e > esyn,KN, with |Aq;| reaching ~ 
0.5 — .75 for q > 30. Hence, for very hard electron injection spectra, with p < 1, the spectral 
index a {— 0.5 — |Aq;|) can even reach negative values. The hardening is visible already at 

e ~ 0.1esyn,KN; 

• For 1 < b^ax ^ bg, the hardening of the synchrotron spectrum combined with the high 
energy break at eg^max leads to the formation a "bump" in the high-energy portion of the 
synchrotron spectrum. For bmax ^ bg, the hardening of the synchrotron spectrum stops at 
esyn,s and the spectrum continues with the same slope as in the Thomson regime but with a 
normalization q times larger; 

• For hard electron injection functions (p < 2) and bmax > bg, the luminosity of the syn- 
chrotron component is larger than luminosity of the IC component, even for g S> 1. This is 

because for a hard injection function most of the power is supplied to electrons with 7 > 7^, 
and the energy losses for these electrons arc dominated by synchrotron radiation. For p = 2, 
the luminosities of the synchrotron and IC spectral peaks are of the same order; 

• When KN effects are important, both the IC and especially the synchrotron component 
can have spectra harder than the hardest spectrum possible in the Thomson limit for fast 
cooling electrons {a = 0.5); 
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• Generically, photon-photon pair production of the Compton gamma-rays on the ambient 
radiation field may be important for •jmax Ikn- However, for a pair production energy 
threshold exceeding e/c,s, the fraction of the 7-ray luminosity converted into pairs is not 
significant, even if the opacity for pair production is large; 

• The continuous energy loss approximation for the evolution of the electron distribution 
appears to work reasonably well, even for ^max ^ Ikn- Use of this approximation can save 
considerable computing time; 

• The KN effects we have discussed can be important in powerful blazars and HMXB, with 
the latter including accreting compacts objects and rotationally powered pulsars. 
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A. APPENDIX 

A.l. Scattering of directed photon beams on isotropically distributed 

relativistic electrons 

The general formula for the distribution in energy and angle of the scattered photons 
per electron per unit time for the photon beams is given by Eq. (14) in Aharonian & Atoyan 
(1981). For e 3> eo and 7 3> 1, this formula takes the form (see Eqs.[20] and [21] in Aharonian 
& Atoyan 1981) 

dNsc{eri.O) 3 f n. 



f Ti 

car / —^f{e,eo,^,0)deo, (Al) 



dedn 167r 

where eo and e are energies of the incident and scattered photons in rrieC^ units, respectively, 
9 is the scattering angle, is the photon number density per energy, 

'''-^^^ = 2(l-cosg)7^(l-(e/7)) ' ^^'^ 

and 

/(e, eo, 7, ^) = 1 + 7^. r - n + I271 w (A3) 

2(1— b0[l — w) bQ[l — wy 

where Bq — 2(1 — cos^)eo7, and w — e/^. 
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A.2. Scattering of isotropically distributed photons on isotropically 

distributed relativistic electrons 



For an isotropic radiation field 



dedn Att de IQir 
where (see Eq.[23] in Aharonian & Atoyan 1981) 



f Tl 

car / —^fiso(e,eo,-f)deo, (A4) 



\ f 2 rcosO„ 



and 

/iso(e,eo,7) = ^ / f{^,^o,1^0)dQ = ^ I f {e,eo, -f,9)d cos 9 , (A6) 

where cos 6min = 1 — 2^^/(6(1 — w)). Integration in Eq. (A6) can be performed analytically 
giving (see Eq.[22] in Aharonian & Atoyan 1981) 

. , , ^ w'^ w 2w'^ 2w b(l — w) 

tisoi^, Cn, 7) = 1 H ; + m , 

^ ^ 2(1 -w) b{l-w) b^{l-wy 26(1 -w)2 b{l-w) w 

(A7) 

where b — 4eo7 (note that fiso is fully equivalent to the term bracketed in Eq.[9] in Jones 
[1968]). 



A. 3. Electron energy losses 



The rate of inverse- Compton energy losses of electrons is 



fiso{e, eo,7)ede 



den . 



l7kc - -^car-, J - 
The inner integral has an analytical solution (see Eq. [46] in Jones 1968) 



J fiso{e, eo,-f)ede 



where 



^6 + 6 + £^ ln(l + 6) - (^^b^ + 66^ + 96 + 4 
and Li2 is the dilogarithm. Hence, 



-2 + 2Li2(-6) 



n\ic 



^^7^ / fKN{b)eoneodeQ 



where 



fKN{b) = 9g{b)/b' 



(A8) 

(A9) 
(AlO) 

(All) 
(A12) 
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Fig. 1. — The function FxN{b) computed for mono-energetic ("mono") and power law (ag = 
0.0 and ao = 0.5) energy distributions of the external photon field. The solid lines show the 
results of the exact calculations while the dashed lines are the analytical approximations. 
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Fig. 2. — bg vs q for: mono-energetic and power law ambient radiation fields. Solid lines - 
exact results, dashed lines - analytical aproximations. 
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Fig. 3. — The relative IC and synchrotron energy losses for an ambient radiation energy 
distribution that is mono-energetic (solid lines) or power law (cto = 0.0 - dotted lines; 
ao — 0.5 - dashed lines). 
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Fig. 4. — Inelasticity, A = e/c(7)/7 fo^ electron Compton scatterings off a mono-energetic 
ambient radiation field. Dotted line - isotropic ambient radiation field. Solid lines - beamed 
ambient radiation field, for scattering angles cos^ = — 1.0, 0.0,and0.6. 
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Fig. 5. — Steady state electron energy distributions for power-law electron injection function, 
Q oc 7~P, and mono-energetic ambient radiation field. Solid lines - exact results; dotted lines 
- results obtained using the continuous energy loss approximation for all Compton scattering. 
The model parameters are: p = 2, g = 30; hmax = 1; 10; 10^; 10^. 
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Fig. 6. — Inverse- Compton plus synchrotron spectra of steady sources, with model param- 
eters p = 2; g = 30; bmax = 1, 10, 10^, 10^, 10^; eo = 10"^; B = IGauss. (a) Upper panel - 
mono-energetic ambient radiation field, (b) Lower panel - power-law ambient radiation field 
with ao = 0.5. Solid lines - exact calculations. Dotted lines - calculations using continuous 
energy loss approximation. Dashed lines - Compton spectra computed using the continuous 
energy loss approximation and the delta-function approximation. The dot-dashed line in the 
lower panel is the asymptotic power law (a = —0.5) for the IC spectrum at e > eic,s given 
by Eq. (31). To show convergence to this spectrum for increasing bmax, the lower panel also 
shows the IC spectra obtained for bmax = 10^, 10^, and 10''. 
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Fig. 9. — Inverse-Compton spectra for a beamed ambient radiation field. Tlie model pa- 
rameters: q — 30, bmax = 10, 10"^; eo = 10"*^; cos6' = —1.0,0.0,0.6. For comparison, the IC 
spectra for an isotropic ambient radiation field are also shown (dotted lines). 




Fig. 10. — Schematic illustration showing the location of the high-energy break in the Comp- 
ton spectra as a function of q. (For details, see text.) 



